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ABSTRACT 

The CG4 and SalOl regions together cover a region of ^0.5 square degree in the vicinity of a 
"cometary globule" that is part of the Gum Nebula. There are seven previously identified young 
stars in this region; we have searched for new young stars using mid- and far-infrared data (3.6 
to 70 microns) from the Spitzer Space Telescope, combined with ground-based optical data and 
near-infrared data from the Two-Micron All-Sky Survey (2MASS). We find infrared excesses in 
all 6 of the previously identified young stars in our maps, and we identify 16 more candidate 
young stars based on apparent infrared excesses. Most (73%) of the new young stars are Class II 
objects. There is a tighter grouping of young stars and young star candidates in the SalOl region, 
in contrast to the CG4 region, where there are fewer young stars and young star candidates, and 
they are more dispersed. Few likely young objects are found in the "fingers" of the dust being 
disturbed by the ionization front from the heart of the Gum Nebula. 

Subject headings: stars: formation - stars: circumstellar matter ~ stars: pre-main sequence - 
infrared: stars 



1. Introduction 

Hawarden and Brand (1976) identified "several elongated, comet-like objects" in the Gum Nebula. 
These objects have dense, dark, dusty heads and long, faint tails, which are generally pointing away from 
the center of the Vela 0B2 association. More such "Cometary Globules" (CGs) were subsequently identified 
in the Gum Nebula (e.g., Sandqvist 1976, Reipurth 1983), but similar structures had also been identified 
elsewhere (e.g., Orion, Rosette Nebula, etc.) in the context of Bok Globules (e.g., Bok & Reilly 1947) and 
"elephant trunks" (e.g., Osterbrock 1957). These objects are all thought to be related in the following sense 
- certain regions of the molecular cloud are dense enough to persist when the stellar winds and ionizing 
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radiation from OB stars powering an H II region move over them, initially forming elephant trunks and 
then eventually cometary globules. These structures often also have bright rims, thought to originate from 
the OB stars' winds and radiation, and are often actively forming stars (e.g., Reipurth 1983), most likely 
triggered by the interaction with winds and radiation from the OB stars (e.g., Haikala et al. 2010). 

Cometary Globule 4 (CG4) in the Gum Nebula has a striking appearance (see Figures [THl]) , where the 
combination of the original dust distribution plus ablation from the OB winds and radiation has resulted in 
a relatively complicated structure. A serendipitous placement of background galaxy ESO 257- G 019 just 
0.15° to the East of the heart of CG4 adds to the drama of the image (see Figures [TH5]) . About a half a 
degree to the West of the heart of CG4 is another cloud, named SalOl. This region was initially recognized 
by Sandqvist (1977) as an opacity class 5 (on a scale of 6, e.g., fairly dark) dark cloud. This cloud appears 
to have been shadowed, at least partially, by CG4 from the ionization front. (See discussion in Pettersson 
2008.) 

Reipurth & Pettersson (1993) studied the region of CG4+Sal01, finding several Ha emission stars; see 
Table [1] They point out that CG-Hal and 7 are not associated with dusty material, and, as such, may 
have been associated with clumps that have already evaporated, as opposed to the stars still projected onto 
molecular cloud material. They argue on the basis of Ha equivalent widths that these cannot be foreground 
or background dMe stars, but are instead likely young stars, members of the association. They then use 
this to argue that this cloud is most likely a part of the Gum Nebula Complex (as opposed to a foreground 
or background object). We therefore assume as well that the stars associated with CG4+Sal01 are also 
associated with the Gum Nebula. 

Distances to CG4-|-Sal01 and even the Gum Nebula are uncertain, with values between 300 and 500 
pc appearing in the literature (e.g.. Franco 1990). The generally accepted source of strong ultraviolet (UV) 
radiation is 7^Velorum; which is taken to be 360-490 pc away (Pozzo et al. 2000). However, the Gum Nebula 
is elongated along our line of sight, so the distance to different parts of the nebula could be significantly 
different than the distance to 7^Velorum. Vela 0B2 is ^ 425 pc (Pozzo et al. 2000). In the context of this 
paper, we test the extrema of the distance estimates of 300 and 500 pc, though we note that our results are 
not strongly dependent on distance. 

Since the CG4-|-Sal01 region contains some previously identified young stars, it is likely that there are 
more young stars, perhaps lower mass or more embedded than those discovered previously. Kim et al. (2003), 
using a preliminary reduction of some of the optical data used here, discussed some additional candidate 
young stars in this region. Since it is now commonly believed that every low-mass star goes through a period 
of having a circumstellar disk, young stars can be identified via an infrared (IR) excess, assumed to be due 
to a circumstellar disk. A survey in the IR can be used to identify objects having an IR excess and thus 
distinguish candidate young stars from most foreground or background objects, at least those foreground 
or background stars without circumstellar disks. The IR also more easily penetrates the dusty environs of 
star-forming regions, particularly globules such as these cometary globules in the Gum Nebula. 

The Spitzer Space Telescope (Werner et al. 2004) observed the CG4-|-Sal01 region with the Infrared 
Array Camera (IRAC; Fazio et al. 2004) at 3.6, 4.5, 5.8, and 8 fim, and with the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) at 24 and 70 fim. We used these data to search this region 
for additional young stellar object (YSO) candidates. We combined these Spitzer data with data from the 
near-infrared Two-Micron All-Sky Survey (2MASS; Skrutskie et al. 2006) and from ground-based optical 
photometric data that we have obtained, and used the multi-wavelength catalog to evaluate and rank our 
list of Spitzer-selected YSO candidates. 
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The observations and data reduction are described in We select YSO candidates using Spitzer colors 
in f]31 and discuss their overall properties in [21 We include a few words on the serendipitously observed 
galaxy in Sj5l Finally, we summarize our main points in [JS] 

2. Observations, Data Reduction, and Ancillary Data 

In this section, we discuss the IRAC and MIPS data acquisition and reduction. We briefly summarize 
the optical (BVRcIc) data reduction, which will be covered in more detail in Kim et al., in preparation. We 
also discuss merging the photometric data across bands, and with the 2MASS near-IR catalog (JHKs). The 
regions of the sky covered by IRAC, MIPS, and the optical observations are indicated in Figure [T] 

We note for completeness that the four channels of IRAC are 3.6, 4.5, 5.8, and 8 microns, and that the 
three channels of MIPS are 24, 70, and 160 microns. These bands can be referred to equivalently by their 
channel number or wavelength; the bracket notation, e.g., [24], denotes the measurement in magnitudes 
rather than flux density units (e.g., Jy). Further discussion of the bandpasses can be found in, e.g., the 
Instrument Handbooks, available from the Spitzer Science Center (SSC) or the Infrared Science Archive 
(IRSA) Spitzer Heritage Archive (SHA) websites. 

2.1. IRAC Data 

We used the IRAC data for CG4 from program 462, AORKE^ 24250880; for SalOl, we used the IRAC 
data from program 20714, AORKEY 16805888. The IRAC data from program 462 (for CG4) were taken 
on 2007-12-29 with 12 sec high-dynamic-range (HDR) frames, so there are two exposures at each pointing, 
0.6 and 12 sec, with 3 dithers per position, for a total integration time of 36 seconds (on average). The 
IRAC data from program 20714 (for SalOl) were taken on 2006-03-27 with 30 sec HDR frames; there are 
also two exposures per pointing, but deeper, 1.2 and 30 seconds. For this observation, there are two dithers 
per position, for a total integration time of 60 seconds (on average). Because of the different integration 
times, we reduced the SalOl and CG4 observations independently even though they overlap on the sky (in 
Figures [^HSJ the jagged-edged observation is from the observation in program 462, and the smooth-edged 
observation on the right is from the observation in program 20714). 

We note that there are additional IRAC data in this region that we did not use. IRAC data from 
program 202 were of a very small region centered on the head of the globule. We did not include these data 
in an effort to make our survey as uniform as possible over the entire surveyed region. IRAC data from 
program 20714 for CG4 were taken in non-HDR mode, with 30 second exposures; as a result of some very 
bright stars in the field of view, the instrumental effects rendered these data very difficult to work with. Our 
science goals near CG4 can be met with the total integration time from program 462 alone. In addition, the 
data from program 462 cover a larger area (~ 0.8°x ~ 1°) than from program 20714 (~ 0.4°x ~ 0.5°). For 
these reasons, we did not incorporate the CG4 IRAC data from program 20714 in this analysis. 

We started with the corrected basic calibrated data (CBCDs) processed using SSC pipeline version 18.7. 
Because of the very bright stars in the field of view, and because the data from program 462 were taken 



^An AOR is an Astronomical Observation Request, the fundamental unit of Spitzer observing. An AORKEY is the unique 
8-digit identifier for the AOR, which can be used to retrieve these data from the Spitzer Archive. 
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Fig. 1. — Approximate location of optical, IRAC, MIPS coverage, superimposed on a Palomar Observatory 
Sky Survey (POSS) image. Purple is optical, blue is IRAC-1 and -3, green is IRAC-2 and -4, orange is 
MIPS-1, and red is MIPS-2. The approximate locations of CG4 and SalOl are also indicated. The galaxy 
ESO 257- G 019 is located on the left, partially obscured by the edge of the optical coverage. 
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Fig. 2. — The IRAC 3.6 /im (channel 1) mosaic. 
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Fig. 3. — The IRAC 4.5 /im mosaic (channel 2). 
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Fig. 4. — The IRAC 5.8 fim mosaic (channel 3). 




Fig. 5. — The IRAC 8 mosaic (channel 4). 
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with cluster targets, we could not use the pipeline-processed mosaics. Moreover, the artifact correction, 
which is normally done for individual cluster targets separately in the SSC pipeline processing, is much 
improved when using the CBCD files from program 462 all at once. We reprocessed the IRAC data from 
both program 462 (for CG4) and program 20714 (for SalOl), using MOPEX (Makovoz & Marleau 2005) to 
calculate overlap corrections and create mosaics with very much reduced instrumental artifacts compared to 
the pipeline mosaics. The pixel size for our mosaics was the same as the pipeline mosaics, 0.6 areseconds, 
half of the native pixel scale. We created separate mosaics for the long and the short exposures at each 
channel for photometric analysis. For display purposes, we further used MOPEX to combine the two long- 
frame observations into one large mosaic per channel, as seen in Figures The component mosaics were 
properly weighted in terms of signal-to-noise and exposure time. The total area covered by at least one 
IRAC channel (as seen in Figures [211S]) is '-^0.5 square degrees. 

To obtain photometry of sources in this region, we used the APEX-lframe module from MOPEX to 
perform source detection on the resultant long and short mosaics for each observation separately. We took 
those source lists and used the aper.pro routine in IDL to perform aperture photometry on the mosaics 
with an aperture of 3 native pixels (6 resampled pixels), and an annulus of 3-7 native pixels (6-14 resampled 
pixels). The corresponding aperture corrections are, for the four IRAC channels, 1.124, 1.127, 1.143, & 1.234, 
respectively, as listed in the IRAC Instrument Handbook. As a check on the photometry, the educators and 
students associated with this project (see Acknowledgments) used the Aperture Photometry Tool (APT; 
Laher et al. 2011a, b) to confirm by hand the measurements for all the targets of interest. To convert the 
flux densities to magnitudes, we used the zero points as provided in the IRAC Instrument Handbook: 280.9, 
179.7, 115.0, and 64.13 Jy, respectively, for the four channels. (No array-dependent color corrections nor 
regular color corrections were applied.) We took the errors as produced by IDL to be the best possible 
internal error estimates; to compare to flux densities from other sources, we took a flat error estimate of 5% 
added in quadature. 

To obtain one source list per channel per observation, we then merged the short and the long exposures 
for each channel separately, and for each observation independently because of the different exposure times 
as noted above. The crossover point between taking fluxes from the short and long exposures were taken 
from empirical studies of prior star-forming regions, and were magnitudes of 9.5, 9.0, 8.0, & 7.0 for the 
four IRAC channels repsectively. We performed this merging via a strict by-position search, looking for 
the closest match within 1". This maximum radius for matching was determined via experience with other 
star- forming regions (e.g., RebuU et al. 2010). The limiting magnitudes of these flnal source lists are the 
same for both observations, and are [3.6]~17 mag, [4.5]~17 mag, [5.8]'--^15.5 mag, and [8]~14.5 mag. 

2.2. MIPS Data 

There are two MIPS AORs in the CG4 region and two MIPS AORs in the SalOl region, all four of which 
were obtained as part of program 20714 (AORKEYs 16805632, 16807936, 16806144, 16808192) on 2006-05- 
08 or 2006-06-12; see Figures [TJ [51 and[71 The AORs were obtained in MIPS photometry mode, nominally 
centered on 7:31:18.7, -46:57:45 for SalOl and 7:33:48 -46:49:59.9 for CG4. One of each pair of AORs is 
explicitly a MIPS-24 photometry observation, and the other is a MIPS-70 photometry observation. During 
the 70 fj,m observation, the 24 ^m array is still turned on and is still collecting valid data. We combined 
the prime 24 ^m data from the MIPS-24 photometry mode observations with the serendipitous 24 ^m data 
from the MIPS-70 photometry mode observations to obtain larger maps at 24 /im. The original explicitly 
24 /im photometry observations are small photometry-mode maps, with 3 s integration per pointing, but 
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Fig. 6. — The MIPS 24 /j,m mosaic. The 24 /im coverage consists of two pointed photometry-mode small 
maps, plus the 24 /im data serendipitously obtained during 70 fim photometry observations. Extended 
emission and point sources are both apparent. 
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Fig. 7. — The MIPS 70 /zm mosaic. The 70 /xm coverage is two pointed photometry-mode small maps. 
Extended emission and point sources are both apparent. 
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a net integration time of ~210 s over most of the resultant mosaic. The screndipitously obtained data 
averaged ^350 s net integration time; where the dehberate and serendipitous data overlapped, integration 
times can be '^450 s. The 70 fiia photometry observations were also small photometry-mode maps, with 10 
s integration per pointing; the net integration time over most of the mosaics was ^400 s. 

The data for 24 /im, like those for IRAC, were affected by the bright objects, and required additional 
processing beyond what the online pipelines could provide. We started with S18.13 enhanced BCDs (eBCDs) 
from the pipeline. We implemented a self-flat for each AOR separately, as described in the MIPS Instrument 
Handbook, available from the SSC or IRSA SHA websites. For each pair of overlapping 24 micron maps, 
we then ran an overlap correction using the overlap script that comes with MOPEX, and then created one 
mosaic for CG4 and one for SalOl, again using MOPEX. Our mosaics had the same pixel size as the online 
mosaics, 2.45". In order to combine the images into one mosaic for display in Figure El the different overall 
background levels between the two observations (having an origin in the different Zodiacal light levels at the 
times of the two observations) were problematic. The brighter of the two was artificially lowered via median 
subtraction to bring its dynamic range into a similar regime as the fainter; this renders photometry on this 
net mosaic invalid, but the morphology seen in Figure IH] is still valid. The total area covered by the net 24 
fim map is only '^0.3 square degrees, smaller than that of the IRAC map. 

To obtain photometry at 24 /im, we ran APEX-lframe on each of the mosaics (one per observation) and 
performed point-response-function (PRF) fitting photometry using the SSC-provided PRE. Tests using the 
apex_qa module portion of MOPEX suggest that our photometry is well within expected errors. For three 
problematic sources, we used aperture photometry instead of the PRE-fitted photometry, as they provided 
a better fit in apex_qa. We used the signal-to-noise ratio (SNR) value returned by APEX-lframe as the best 
estimate of the internal (statistical) errors, adding a 4% flux density error in quadrature as a best estimate 
of the absolute uncertainty. The limiting magnitude of these observations is [24]~10.5 mag. Note that we 
optimized our data reduction to obtain measurements of the brighter sources and sources superimposed on 
the nebulosity; many sources fainter than this are apparent in the image but not included in our catalog, 
simply because our scientiflc goals are aimed at the brighter objects. For one source of interest below 
(073243.5-464941, which was considered and then rejected as a YSO candidate; see 0, an upper limit was 
obtained at the given position by laying down an aperture as if a source were there, and taking 3 times that 
value for the 3cr limit. To convert the flux densities to magnitudes, we used the zero point as found in the 
MIPS Instrument Handbook, 7.14 Jy. 

At 70 /im, there are viable observations from AORKEYs 16807936 and 16808192. We downloaded data 
processed with pipeline version S18.12. The online pipeline does a very good job of producing mosaics; see 
Figure [3 where there are a handful of point sources and extended emission visible. The online pipeline 
produces both flltered and unfiltered mosaics; the filtering preserves the flux densities of the point sources 
and improves their signal-to-noise, especially for faint sources, but destroys the flux density information for 
the extended emission. The unfiltered mosaics are shown in Figure [71 but we performed photometry on 
the filtered mosaics. The pipeline mosaics have resampled 4" pixels (as opposed to 5.3" native pixels), and 
the two observations together cover about 0.1 square degrees. We used APEX-lframe to do PRF fitting 
on the pipeline filtered mosaics for the point sources, using the SSC-provided PRF. For one problematic 
source, aperture photometry provided a better fiux density estimate. There are only 11 objects with 70 ^m 
detections, and there is a large variation in background levels, so quoting a limiting magnitude is difficult, 
but is very approximately 3 mag. We assumed a conservative, fiat 20% fiux density error. The zero point 
we used again came from the MIPS Instrument Handbook, 0.775 Jy. 

Where there was 70 /im coverage for the sources of interest, we placed an aperture at the expected 
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location of the source, and performed photometry as if there were a source there, taking 3 times that vahie 
as the 3cr hmit that appears in Table [5] below. 



2.3. Optical Data 



The optical data will be discussed further in Kim et al., in preparation, but we summarize the important 
aspects of the data reduction here. 

The BVRcIc photometry of the CG4+Sal01 region were obtained during 2001 March 6, 7, 9, 10, and 
11 using the 2Kx2K CCD at the 0.9m telescope at the Cerro Tololo Inter-American Observatory (CTIO). 
The images have a pixel scale of 0.4"in a 13.6'field of view. 

Bias and twilight sky flat fields were taken at the beginning and at the end of each night. Long and short 
(300 sec and 30 sec) exposures were taken for object fields. During every photometric night we observed 
Landolt (1992) standard stars of two or three fields several times per night for photometric calibration. 

We performed aperture photometry using multiple aperture sizes, and the photometry with highest 
signal-to- noise (S/N) ratio were chosen as final photometry for each star. For the standard stars, we used 
an aperture radius of 17 pixels. We used IRAfH/PHOTCAL routines to solve for the zero point, extinction 
and color terms of the standard star solution. 

For the target stars in CG4 and SalOl, we used a custom IDL photometry pipeline (written by 
W. H. Sherry), which was developed for the CTIO 0.9m telescope. For each target, aperture photome- 
try was performed using multiple size apertures starting from aperture size of 2 pixels to 17 pixels. The 
highest S/N photometry was chosen as the final photometry. 

We used an aperture correction to place our photometry on the same photometry system as our standard 
stars. The point-spread function (PSF) of the CTIO 0.9m telescope varies noticeably as a function of location 
on the CCD. This is insignificant for large apertures, but for aperture sizes of 2-3 pix, the difference can be a 
few percent. We accounted for the spatial dependence of the aperture correction in each image by fitting the 
aperture corrections for stars with photometric errors better than 0.02 mag with a quadratic function. The 
uncertainty of the aperture correction is about 0.01 magnitude. For each aperture, the uncertainty on the 
instrumental magnitude is calculated including the uncertainty on the aperture correction for each aperture. 

We used imwcs program written by D. Minl(f| (Mink 1997) to determine the astrometric solution of each 
image. It fits the pixel coordinates of our targets to the known positions of USNO A2.0 stars located in each 
field. We then measured positions to an accuracy of ^0'.'3 relative to the USNO A2.0 reference frame. 

For each pointing, we used a 1" matching radius to match sources within the optical BVRcIc filters. 
The coordinates are averaged over different filters, and the typical average uncertainty is 0'.'2~0'.'3. To find 
duplicates between adjacent pointings, we again looked for positional matches within 1", and then took a 
weighted mean of the available photometry. 

Basically the entire IRAC map was covered west of 113.8 degrees RA (07:35:12); see Figure[TJ We have 
found no YSO candidates east of this. The completeness limits over the field are as follows: B ^ 19 mag. 



^IRAF is distributed by the National Optieal Astronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 

^Documentation and source code are available at |http:/ /tdc- www. harvard. edu/software/wcstools/| 
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V ^ 18 mag, Rc ^ 17.5 mag, Ic 17 mag. However, we note that there are fewer sources per projected 
area detected (e.g., effectively shallower limits) in the regions where there is molecular cloud material. The 
zero-points we used for conversion of the magnitudes to flux densities (for inclusion in the spectral energy 
distributions in Figures [HHH]) were, respectively, 4000.87, 3597.28, 2746.63, and 2432.84 Jy. 

2.4. Bandmerging 

In summary, to bandmerge the available data, we first merged the photometry from all four IRAC 
channels together with near-IR 2MASS data within each observation (CG4 and SalOl), and then merged 
together those source lists from each observation. We next included the MIPS data, and then the optical 
data. We now discuss each of these steps in more detail. At the end of this section, we discuss some aggregate 
statistics of the bandmerged catalog. 

To merge the photometry from all four IRAC channels together, we started with a source list from 
2MASS. This source list includes JHKg photometry and limits, with high-quality astrometry. We merged 
this source list by position to the IRAC-1 source list, using a search radius of 1", a value empirically 
determined via experience with other star- forming regions (e.g., RebuU et al. 2010). Objects appearing in 
the IRAC-1 list but not the JHKg list were retained as new potential sources. The master catalog was then 
merged, in succession, to IRAC-2, 3, and 4, again using a matching radius of 1". 

Because we are primarily interested in objects detected by Spitzer, we dropped any objects not having 
flux densities in at least one Spitzer band (e.g., objects off the edge of the Spitzer maps, having measurements 
only in 2MASS). Because the source detection algorithm we used can be fooled by instrumental artifacts, 
we also explicitly dropped objects seen only in one IRAC band as likely artifacts. 

For the sources of interest later in the paper, most have counterparts in all three bands in 2MASS by 
this point in the merging. However, for two sources, this matching failed, at least in part. For one source, 
073121.8-465745 (=CG-Ha4), the automatic merging found a counterpart with a measured J magnitude, 
but the HKg measurements are flagged with a photometric quality (ph_qual) flag of 'E', denoting that 
the goodness of fit quality was very poor, or that the photometry fit did not converge, or that there were 
insufficient individual data frames for the measurement. Since there was a good measurement at J, we 
assumed that there were sufficient frames at HKg, and that something had happened to the fit. We took 
the values as reported in the catalog, assumed a large error bar, and used these values in the table and 
plots below; as will be seen in the spectral energy distribution (SED) for this object (Fig. [151 top center), 
these values are probably close to what is most appropriate for this object. For source 073425.3-465409, the 
automatic merging fails, most likely because the 2MASS counterpart is slightly extended, and it may be 
slightly extended at 3.6 /zm as well. The nearest match in the catalog is ~ 2" away, very large by comparison 
to other source matches here, but manual inspection strongly suggests that this is the appropriate counterpart 
to the source seen at Spitzer bands. See Appendix lA. 211 for more on this interesting source. 

We then compared the source lists from the separate observations, CG4 and SalOl, again using a 
matching radius of 1". For objects detected in more than one mosaic, we took a weighted average of the flux 
density at the corresponding band. 

Next, we merged the 70 fim source list to the 24 /im source list. The 70 fim point-spread function is 
large compared to the positional accuracy needed, and astrophysically, each 70 /im source ought to have 
a counterpart at 24 /im, given the sensitivity of these observations. We individually verified that each of 
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the 70 /iiTL point sources had a 24 /^m counterpart, and then merged these two source hsts. To successfully 
have the computer match the sources that were clearly matches by eye, a positional accuracy of 2.5" was 
required, consistent with our experience in other star- forming regions (e.g., Rebull et al. 2010). Since the 
two MIPS observations do not overlap with each other, no explicit merging of the MIPS source lists from 
the two observations was required beyond simple concatenation. 

To combine the merged MIPS source list into the merged 2MASS-I-IRAC catalog, wc used a positional 
source match radius of 2", again determined via experience with other star- forming regions (e.g., Rebull 
et al. 2010). The MOPEX source detection algorithm can be fooled by structure in the nebulosity in the 
image, and by inspection, this was the case for these data. To weed out these false 'sources', we then dropped 
objects from the catalog that were detections only at 24 fim and no other bands. Finally, to merge the J 
through 70 iim catalog to the optical (BVRcIc) catalog, we looked for nearest neighbors within 1". 

After this entire process, there are ^21,000 sources with IRAC-1 (3.6 ^m) or IRAC-2 (4.5 fim) detections, 
~9000 sources with IRAC-3 (5.8 iim) detections, and ^4000 sources with IRAC-4 (8 /im) detections. About 
3000 (^^15%) of the IRAC sources have viable data at all 4 IRAC bands, nearly all of which have counterparts 
in 2MASS. The optical data do not cover the entire IRAC map, but about half of the 4-band IRAC detections 
have counterparts in the optical catalog. There are only '-^500 sources at 24 jim in our catalog and just 11 
sources at 70 /im; note that the MIPS-24 map covers a much smaller area than the IRAC maps, and the 
MIPS-70 map is smaller still (see Figures [l][7]) . Ten of the 11 MIPS-70 sources have counterparts at all four 
IRAC bands; the one that docs not is saturated at 2 of the IRAC bands. About 200 sources have all four 
IRAC bands plus MIPS-24. 

3. Selection of YSO candidates w^ith infrared excess 

With our new multi- wavelength view of the CG4-|-Sal01 region, we can begin to look for young stars. We 
focus on finding sources having an infrared excess characteristic of YSOs surrounded by a dusty disk. There 
is no single Spitzer color selection criterion (or set of criteria) that is 100% reliable in separating members 
from non- member contaminants. Many have been considered in the literature (e.g., Allen et al. 2004, Rebull 
et al. 2007, Harvey et al. 2007, Gutermuth et al. 2008, 2009, Rebull et al. 2010, Rebull et al. 2011). Some 
make use of just MIPS bands, some make use of just IRAC bands, most use a series of color criteria, and 
where possible, they make use of (sometimes substantial) ancillary data. In our case of the CG4-|-Sal01 
region, we have some ancillary data, but the bulk of the data are IRAC+2MASS data. In this case, the 
best choice for selecting YSO candidates is the approach developed by Gutermuth et al. (2008, 2009) and 
adapted by Guieu et al. (2009, 2010) for the case in which no extinction map is available. This selection 
method starts from the set of objects detected at all four IRAC bands and uses 2MASS and MIPS data 
where possible. It implements a series of color cuts to attempt to remove contaminants such as background 
galaxies and knots of nebulosity. 

When we impose these IRAC-based color cuts, we find 25 potential YSO candidates. Wc then inspected 
each of these in all available images and color-color and color-magnitude diagrams. On the basis of this 
inspection, we dropped three of the 25 potential YSO candidates off of our list, though additional data will 
be needed to be sure that these objects are extragalactic. Two of those dropped objects (073542.2-470126 
and 073548.5-470727) have no available data other than IRAC, their SEDs are very flat, and they are located 
near the edges of our images, far from other YSOs and nebulosity. We suspect that these arc extragalactic 
contaminants. The third object, 073243.5-464941, is returned by the IRAC selection as having a small excess 
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Fig. 8. — [3.6] — [4.5] vs. [5.8] — [8] color-color diagram for CG4-|-Sal01. Small dots are objects in the catalog; 
larger dots are objects identified as contaminants, and large red diamonds highlight our YSO candidates. 
All of our Spitzcr-sclccted YSO candidates have colors in this diagram consistent with known YSOs, but 
many contaminants do too. 
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at 5.8 and 8 /im. It is seen at 2MASS and IRAC bands, but is undeteeted at 24 fim to a fairly stringent limit 
(971 ^Jy, or 9.67 mag). If it has an excess at 24 /im, it is a very small excess. The 70 /im data do not cover 
this object, so there are no constraints (not even limits) at 70 fim. Moreover, it is a relatively faint source 
next to a very bright source, located far from any nebulosity. The wings of the bright source arc likely to 
adversely affect the photometric accuracy of the measurements associated with this object. Because of this 
uncertainty and the very low excess as measured, we have dropped this object from our YSO candidate list 
as a likely foreground or background star. 

The remaining 22 YSO candidates that pass the color cuts are shown in Figure[5] In this Figure, objects 
with zero color are likely foreground or background stars (photospheres without disks), though some could 
be young stars that have already shed their disks. Young stars with circumstellar disks are generally red 
in both IRAC colors, but contaminants such as galaxies also may have these colors. All of the 22 objects 
highlighted in this Figure have IRAC colors consistent with young stars with disks. 

The Gutermuth et al. (2008, 2009) selection criteria have provisions for adding stars to the list of 
candidate YSOs based on properties at other bands, such as MIPS bands. We now investigate the properties 
of objects in our catalogs at the MIPS bands to see if we should add additional objects to our list of YSO 
candidates. In summary of the rest of this section, while we find some interesting objects, in the end, we do 
not add any more YSO candidates to our list. 

Young stars having inner disk holes and thus excesses at only the longest bands can be revealed via 
comparison of the 24 /im measurement to a shorter band, such as Kg or [3.6]. If the data are available, 
one should use [3.6] vs. [3.6] — [24] rather than Kg vs. Kg — [24]. There is an intrinsic spread in Kg — [24] 
photosphcric colors that is not present in [3.6] — [24] because late type stars are not colorless at Kg — [24] 
(Gautier et al. 2007). The effects of reddening are stronger at Kg than at 3.6 fj,m. And, if 2MASS is the 
only source of Kg, even short 3.6 /im integrations can reach fainter sources than 2MASS does. In our case 
of CG4+Sal01, the IRAC coverage is larger than the MIPS coverage, and so we use [3.6] vs. [3.6] — [24] to 
look for any objects with an excess starting at 24 fim. 

Figure [9] shows this [3.6] vs. [3.6] — [24] diagram, with the same notation as the prior figure. Ordinary 
stellar photospheres (likely foreground or background stars) have [3.6] — [24] ~0, and galaxies make up the 
large, elongated source concentration near [3.6] — [24]^6, [3.6]^16. Objects not in this region, e.g., the brighter 
and/or redder objects, are less likely to be part of the Galactic or extragalactic backgrounds, and more likely 
to be YSOs with a 24 /xm excess. Most of the IRAC-identified YSO candidates are indeed in the region of 
this diagram occupied by other known YSOs (see, e.g., Rebull et al. 2010, 2011, Guieu et al. 2009, 2010). 
There is one (073049.1-470209) that is among the reddest objects in this diagram, near [3.6] — [24]^10; see 
Appendix ?? for more on this specific object. Most of the objects already ruled out as YSO candidates based 
on their IRAC properties are in the extragalactic concentration of sources. The objects with [3.6] — [24] ^^0 
do not have apparent excesses, but there are eight additional objects with [3.6] — [24]>1 and [3.6]<14.5 that 
seem to have the right placement in this diagram to be YSO candidates. We investigated each of these 
candidates, and none had evidence based on SED shape, significance of excess, or appearance in the images 
compelling enough to have us add them to our list of YSO candidates. The apparent small excesses just at 
24 fim are most likely due to source confusion at the lower resolution 24 /xm band, with either a background 
source or a low-mass companion. The most compelling one based on the SED is 073355.0-464838, but the 
source seems to be confused with a nearby source that emerges at 8 /im, and we strongly suspect that the 
24 /im flux instead corresponds to the object appearing at 8 /im, rather than the point source seen at 8 /im 
and shorter bands. We do not add this source to our list of YSO candidates at this time. Higher spatial 
resolution 24 /im observations would be required to resolve this issue. 
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Fig. 9. — [3.6] vs. [3.6] — [24] color-magnitude diagram for CG4+Sal01. Small dots are objects in the catalog; 
larger dots are objects identified as contaminants, and large red diamonds highlight our IRAC-selected YSO 
candidates. All of the IRAC-selected YSO candidates have colors in this diagram consistent with known 
YSOs. 
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Fig. 10. — [24] vs. [24] — [70] color-magnitude diagram for CG4+Sal01. As in prior plots, dots are objects 
identified as contaminants, and large red diamonds highlight our IRAC-selected YSO candidates. A red 
square highlights a very bright likely background star; see the text. The arrows indicate the positions of the 
YSO candidates for which we could obtain upper limits at 70 fim. The IRAC-selected YSO candidates have 
colors in this diagram consistent with known YSOs. 
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Some of the brightest stars in the CG4+Sal01 region are saturated in at least the first two IRAC 
bands, so neither of the YSO search mechanisms we have used thus far would find them. However, these 
sources could be YSOs, and they are not all saturated in MIPS bands. As our last attempt to search for 
YSO candidates with infrared excesses. Figure [TU] shows the [24] vs. [24] — [70] diagram for our region. This 
diagram for this region is sparse, but a better-populated diagram (see, e.g., Rebull et al. 2010) would basically 
resemble Figure [9l with photospheres being bright and colorless (having [24] — [70]~0), and galaxies red and 
faint. In Figure [101 the few points that are available include things that are galaxies (things that have been 
ruled out as contaminants based on IRAC) which are in the extragalactic part of the diagram here as well, 
six detections that have colors consistent with YSOs, many limits for our YSO candidates, and one thing 
(square in Figure [TO)) that is too bright to be a likely galaxy, and does not appear to be particularly red in 
this diagram. This very bright object is 073339.7-464839, and its SED suggests at first glance that it might 
be a YSO with a small excess just at 24 and 70 /im. It is detected at JHKs. [5.8], and [8], and the Ks-, [5.8], 
and [8] measurements are all consistent with one Rayleigh- Jeans tail, and the [24] and [70] measurements are 
offset on a different, redder Rayleigh- Jeans tail, as might be consistent with a small thermally-heated dust 
disk. However, it is a very bright object, and saturated at the shortest two IRAC bands; the photometry at 
Ji/A's[5.8][8] may also be compromised beyond what our formal errors suggest. It is matched in Simbad to 
IRAS 07321-4642. We suspect that this is a background asymptotic giant branch (AGB) star, or another 
sort of bright background giant, and not a legitimate YSO candidate. We do not include it on our YSO 
candidate list. 

We move ahead from here with the 22 IRAC-selected YSO candidates, and now investigate their multi- 
band properties. 

4. Properties of selected YSO candidates 

4.1. Optical properties 

Optical data can greatly aid in confirming or refuting YSO candidacy because they provide constraints 
on the Wien side of the SED. In Guieu et al. (2010), most of the IRAC-selected candidates in IC 2118 
proved to be too faint, most vividly in the optical, to be likely cluster members. Just ten of our candidate 
YSOs have optical data available, and they appear in Figure [TTJ The objects with optical data that have 
already been ruled out as YSOs based on their IRAC properties are all well below the Siess et al. (2000) 30 
Myr isochrone scaled to 500 pc. One YSO candidate object appears below the 30 Myr isochrone; 073337.6- 
464246 is within the distribution of clear non-YSO points. We do not remove this object from our list, since 
a variety of reasons (such as scattered light) could result in a YSO appearing below the 30 Myr isochone; for 
more discussion of this object, please see Appendix lA. 191 As noted above, the distance to this association is 
uncertain; if the isochrones are instead scaled to 300 pc, then one other object (073121.8-465745, a previously 
identified YSO) appears to be just below the 30 Myr isochone instead of just above it. 

Deeper optical data are desirable in order to obtain magnitude estimates for the remaining YSO candi- 
dates. 
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Fig. 11. — V vs. {V — Ic) color-magnitude diagram for CG4+Sal01. Small dots are objects in the catalog, 
larger dots are objects identified as contaminants, and large red diamonds highlight our YSO candidates. 
Isochrones given are models from Siess et al. (2000) at 1, 10, and 30 Myr, scaled to 500 pc, where we have 
tuned the color-effective temperature relation such that the 100 Myr isochrone matches that of the Pleiades 
single-star sequence (Stauffer et al. 2007, Jeffries et al. 2007). A reddening vector is also indicated. All of 
the YSO candidates shown here, except for one, are in the region occupied by young stars at 300-500 pc. 
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Fig. 12. — J — H vs. H — Ks diagram for the sample, with the same notation as carher figures. The main 
sequence is indicated by a sohd hue. Most of the YSO candidates have an infrared excess starting at H-band 
with moderate reddening. 
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4.2. Near-IR properties 

Near-IR data can also aid in confirming or refuting YSO candidacy. Since we do not liave spectral types 
for most of our sources, it is difficult to estimate the degree of reddening. Figure W2\ shows J — H vs. H — Kg 
for the sample. This plot suggests that most of our YSO candidates have an infrared excess with a moderate 
degree of reddening. 

4.3. B-band properties 

Young stars that are actively accreting from their circumstellar disk can have excess ultraviolet emission 
at [/ or i? bands, or even longer bands during periods of intense accretion. However, these bands are also 
most sensitive to reddening. Figure [13] shows B — V ys. — Ic for the sample, with the main sequence 
indicated as a solid line. Stars with a i?-band excess and relatively small values of Ay would be blue, e.g., 
below the line in this figure. Objects above the line on the upper right of this figure are pushed into that 
location by high Ay- At least 4 of the YSO candidates have a B excess, most likely from mass accretion; 
four more appear to have been pushed from the region of clear B excess by high Ay ■ Similar results for the 
same objects are obtained from the B — V ys. V — Ic plot. The individual objects are listed in the Appendix. 

4.4. Spectral Energy Distributions 

Coordinates and our measured magnitudes between B and 70 /xm for our 22 YSO candidates appear 
in Table [H Six of them (27%) are rediscoveries of the previously known YSOs in this region from Table [1] 
Figures [14]-[11 arc the SEDs for the YSO candidates. 

To guide the eye, we wished to add reddened stellar models to the SEDs, but spectral types are only 
known for the 6 previously known YSOs. In order to provide a reference, for the remaining objects, we 
assumed a spectral type of MO. For each object, a reddened model is shown, selected from the Kurucz- 
Lejeune model grid (Lejeune et al. 1997, 1998) and normalized to band where possible (and to the closest 
band otherwise). Note that this is not meant to be a robust fit to the object, but rather a representative stellar 
SED to guide the eye such that the infrared excesses are immediately apparent. In some cases, ultraviolet 
excesses may also be present. Additional spectroscopic observations are needed to better constrain these fits. 

In the spirit ofWilkinget al. (2001), we define the near- to mid-IR slope of the SED, a = dlog AFA/dlog A, 
where a > 0.3 for a Class I, 0.3 to —0.3 for a flat-spectrum source, —0.3 to —1.6 for a Class II, and < —1.6 
for a Class III. For each of the YSO candidate objects in our sample, we performed a simple ordinary least 
squares linear fit to all available photometry (just detections, not including upper or lower limits) between 
2 and 24 fim, inclusive. Note that errors on the infrared points are so small as to not affect the fitted 
SED slope. The precise definition of a can vary, resulting in different classifications for certain objects. 
Classification via this method is provided specifically to enable comparison within this paper via internally 
consistent means. Note that the formal classification puts no lower limit on the colors of Class III objects 
(thereby including those with SEDs resembling bare stellar photospheres, and allowing for other criteria to 
define youth). By searching for IR excesses, we are incomplete in our sample of Class III objects. The classes 
for the YSO (previously known and candidate) sample appear in Table [S] Out of the 22 stars, 16 (73%) are 
Class II. 

Based on the SEDs and location in several color-color and color- magnitude diagrams, we have ranked 
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Fig. 13. — B — V vs. Rq — Ic for the sample, with the same notation as earher figures. The main sequence is 
indicated by a solid line. Stars with a _B-band excess and relatively small values of Ay would be blue, e.g., 
below the line in this figure. Objects above the line on the upper right of this figure are pushed into that 
location by high Ay. At least 4 and probably 8 of the YSO candidates have a B excess, most likely from 
mass accretion. 



Table 2. Multiband measurements of Spitzer-identified YSO candidates in the CG4+Sal01 region 
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^Values appear in the 2MASS Point Source Catalog with a photometric quality (ph_qual) flag of 'E', denoting that the goodness of fit quality was very poor, or that the photometry fit did not converge. We took the values 
as reported, with a large error bar, as additional constraints on the source. | 

^2MASS JHKs photometry comes from the 2MASS extended source catalog, not the point source catalog; see discussion in sections |2 . 4| and | A . 2 1 1 tO 
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Table 3. Final list of YSO candidates in the CG4+Sal01 region 
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^This grade is meant to indicate rough confidence in the liklihood that the given YSO candidate is a legitimate YSO. Previously identified 
YSOs are given a grade of 'A+\ our highest quality YSO candidates are grade 'A', our mid-grade YSO candidates are grade 'B', and our 
lowest-confidence YSO candidates are grade 'C'. For discussion of individual objects (and an explanation of why each object has that 
grade), please see the Appendix. 
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Fig. 14. — Spectral Energy Distributions (SEDs) for the Spitzer-selected YSOs presented here. + symbols 
are optical data, diamonds are 2MASS (NIR) data, circles are IRAC data, and squares arc MIPS data. If 
a previous source ID exists in Reipurth & Pettersson (1993), it is in the upper right of the plot, along with 
their spectral type (if it exists). Similarly, if there are photometric data from Reipurth & Pettersson (1993), 
they arc marked with light grey dots. The error bars (most frequently far smaller than the size of the symbol) 
are indicated at the center of the symbol. Reddened models from the Kurucz-Lejeune model grid (Lejeune 
et al. 1997, 1998) are shown for reference; see the text. Units of XFx as presented are erg cm~^, and A 
is in microns. 




Fig. 15.— 



SEDs, continued. Notation as in previous figure. 




Fig. 16. — SEDs, continued. Notation as in previous figure. 
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Fig. 17. — Reverse greyscale mosaic of 8 //m data with the locations of the YSOs (previously known and 
new candidates) indicated, color-coded by YSO quality grade. Previously known YSOs are indicated by a 
green dot with a "Y", high-quality (grade 'A') YSOs are indicated by a red dot and an "A", mid-quality 
(grade 'B') YSOs are indicated by a blue dot and a "B", and low-quality (grade 'C') YSOs are indicated by 
a yellow dot and a "C." 



the YSO candidates loosely into three bins: high likelihood of being YSOs (grade A), mid-grade quality 
(grade B), and relatively low likelihood of being YSOs (grade C). This grade also appears in Table [31 Most 
of them (16) are in the grade A bin, which includes the 6 previously identified ones. Comments on individual 
objects (including justifications for the grades that were given) appear in the Appendix. 

Figure [17] shows the 8 /.m mosaic with the positions of the YSO candidates overlaid, color-coded by 
YSO quality. Most of the grade A and B objects are clustered near the previously-known YSOs. The SalOl 
region has a relatively tight clumping of most (16) of the YSOs, with a median nearest neighbor distance 
of ^62". The CG4 region has 6 YSOs, much less tightly clumped, with a median nearest neighbor distance 
nearly 5 times larger, ^301". Clustering is also very commonly found among young stars, so especially in the 
case of the Sal 01 association, the fact that they are clustered also bolsters the case that they are legitimate 
YSOs. 



Because the distance to this association is uncertain, we looked at whether the relative placement of 
stars in the optical color-magnitude diagram could be used to constrain the distance to the stars. Figure [T8l 
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Fig. 18. — Optical My vs. V ~ Ic color-magnitude diagram. The Siess et al. (2000) isochrones are included 
(1, 10, and 30 Myr), but shifted to absolute My. The black stars are our YSO candidates, assuming a 
distance of 500 pc, and the red stars are our YSO candidates, assuming a distance of 300 pc. The grey x 
symbols are Taurus YSOs (from Rebull et al. 2010 and Giidel et al. 2007 and references therein), taken to be 
at 140 pc; Torres et al. (2007, 2009). The Taurus distribution is broad and there are many fewer CG4+Sal01 
stars, but this distribution weakly suggests that CG4+Sal01 is farther rather than closer (see text). 
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presents the My vs. V — I CMD, comparing CG4+Sal01 stars to YSOs in Taurus (with data from RebuU 
et al. 2010, Giidel ct al. 2007, and references therein). Based on morphological grounds (e.g., the degree to 
which the YSOs are still embedded in their natal gas), we expect that the CG4+Sal01 stars might be slightly 
younger than the often more physically dispersed Taurus stars. On the other hand, based on the ratio of 
Class I to Class II sources, the CG4+Sal01 objects might be slightly older than Taurus. In Figure [TSl there 
are not many CG4+Sal01 objects, and the distribution is broad, but assuming a distance of 500 pc, then 
CG4+Sal01 appears to be quite comparable in age to Taurus at '^3 Myr. Assuming a distance of 300 pc, 
the CG4+Sal01 stars are on the whole older than the Taurus stars. Figure [TBI thus weakly suggests that 
CG4+Sal01 is farther rather than closer. 



5. The Galaxy 

The galaxy in our field, ESO 257- G 019, has not been the target of any case studies before. The galaxy 
is classified as type 'SB(s)cd?' in NED, the NASA Extragalactic Database. We measure a major axis length 
of 3f3 or 34.7 kpc at the assumed distance of 36.1 Mpc at the 3a level above the noise in our channel 1 
image. We also measured a minor axis length of 0^8 or 8.4 kpc in the same image for this highly inclined 
galaxy. The surface brightness distribution is close to exponential within the central 3 kpc. The galaxy has 
a clumpy structure at larger distances from its center, with the most prominent clumps appearing about 26" 
and 73" (along the major axis) to the northwest and 26" to the southeast (just outside of the plane of the 
galaxy). We also measured the IRAC [3.6] — [4.5] and [3.6] — [8.0] colors by using all the pixels above the lOtr 
level. The [3.6] — [4.5] color is 0.0 and the [3.6] — [8.0] color is 1.9, both close to the typical values for late-type 
galaxies, as determined by Pahre et al. (2004), and consistent with the galaxy classification as given in NED. 
ESO 257- G 019 appears to be a fairly isolated galaxy, with no nearby companions within our mapped areas. 

6. Conclusions 

We used Spitzer Space Telescope data from the Spitzer Heritage Archive to search for new candidate 
young stars in the CG4+Sal01 region of the Gum Nebula. This region appears to be actively forming young 
stars, perhaps as a result of the passage of an ionization front from the stars powering the Gum Nebula 
(Reipurth & Pettersson 1993). We rediscovered all six of the previously identified young stars in our maps 
as having excesses at Spitzer bands. We have also discovered 16 entirely new young star candidates based 
on their Spitzer properties. We used optical ground-based data and near-infrared data from 2MASS to help 
constrain the SEDs of these new young star candidates. We have sorted the new young star candidates 
into grades of confidence that they are, in fact, legitimate new young stars. We find 16 high confidence 
(grade "A") objects, including the 6 previously identified YSOs, 2 mid-grade confidence (grade "B"), and 4 
low-grade confidence (grade "C" ) young star candidates. For all of the new young star candidates, though, 
additional data will be needed, such as optical photometry where it is missing, and optical spectroscopy to 
obtain spectral types (and rule out extragalactic contamination). Most of the new objects are clustered in 
the SalOl region, and most are SED Class II. 

This work was performed as part of the NASA/IPAC Teacher Archive Research Program (NITARP; 
|http : / / nitarp .ipac.caltech.edu), class of 2010. We acknowledge here all of the students and other folks who 
contributed their time and energy to this work and the related poster papers presented at the January 2011 
American Astronomical Society (AAS) meeting in Seattle, WA. They include: With V. Hoette: C. Gartner, 
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J. VanDerMolen, L. Gamble, L. Matsche, A. McCartney, M. Doering, S. Brown, R. Wernis, J. Wirth, M. 
Bcrthoud. With C. Johnson: R. Crump, N. Kilhngstad, T. McCanna, S. Caruso, A. Laorr, K., Mork, E. 
Stembergs, E. Wigley. With C. Mallory: N. Mahmud. 

We thank J. R. Stauffer for helpful comments on the manuscript. 
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A. Comments on individual objects 



A.l. 



073049.1-470209 



This object, 073049.1-470209, the westernmost in the set, is also the reddest source in the set. It is 
the YSO candidate farthest in the upper right of the IRAC color-color diagram (Fig. [H) and is the reddest 
and faintest YSO in [3.6] vs. [3.6]-[24] (Fig.!}. Based on its SED (Fig. [HI), it is not surprising that there 
are no optical or even 2MASS counterparts; it very steeply falls from 24 /im back to 3.6 ^m, and it seems 
that very deep integrations would be needed to obtain measurements of this object at shorter wavelength 
bands. It is not detected at 70 fim, but it is in a high background region, so the 70 /zm limit is not very 
constraining. This type of SED can be found in very young, very embedded, very low mass objects, but also 
in extragalactic contaminants. Because of its proximity to the rest of the SalOl objects, it remains in the 
list, but we have given it the lowest-quality grade (C). Its steep SED means that it is one of three Class I 
objects in our list. Additional follow-up data are needed to determine the true nature of this source. 



This object, 073049.8-465806, is relatively faint at 3.6 /im, at least in comparison to the rest of the 
YSO candidates here. It appears on the top edge of the clump of likely extragalactic sources in the [3.6] vs. 
[3.6]-[24] plot (Fig. HI). It is only detected at in 2MASS, with limits at J and H. Based on the SED, it 
looks quite extinguished, e.g., there seems to be high Ay in the direction of this source. It is located in close 
(projected) proximity to other YSOs and YSO candidates; YSOs arc also frequently found near other YSOs. 
It is not detected at 70 fim, but the limit is shallow enough that it does not provide a strong constraint 
on the SED. It has an SED characteristic of YSOs, and that plus the apparently high Ay plus its location 
in the cloud close to other YSOs has yielded a "B" grade. The steep rise of the SED from 4.5 to 24 /xm 
influences the classification slope fitting such that it is one of three Class I objects in our list. Additional 
follow-up data are needed. 



Object 073053.6-465742 has 2MASS data, but no available optical data, and just a limit at 70 fim. The 
SED (Fig. [T4| suggests that there is high Ay towards this source. Given our overly simple modelling, there 
seems to be a small but significant excess at 8 /im (better modelling is required to confirm this), and there 
is an apparently large excess at 24 /im. Whenever an apparent excess is seen only at 24 /im, because the 
MIPS-24 camera has lower spatial resolution than the shorter bands, there is a risk that the fiux density 
measured at 24 /im is contaminated by source confusion, most likely with a nearby background source, but it 
could also be with a low mass companion to the young star candidate. There is only a limit at 70 fim to help 
constrain the SED, though it is located in close (projected) proximity to other YSOs and YSO candidates. 
In this case, the fact that the excess appears to be significant at at least 2 bands, its proximity to other 
nearby SalOl sources, and the likelihood that there is high Ay towards this source suggests that it may be 
a legitimate YSO, so we have placed it in the highest quality source bin (grade "A"), and it is a Class II. 
Spectroscopy of this target is required to determine whether or not it is a YSO. 



A.2. 



073049.8-465806 



A.3. 



073053.6-465742 
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A.4. 073057.5-465611=CG-Ha2 

CG-Ha2 was identified in Reipurth & Pettersson (1993) as a YSO. It has counterparts at all bands we 
considered here, including 70 fim. In the optical (Fig. [TT|) . it is in the locus of young stars above the 30 Myr 
isochrone, and it appears there to be the lowest-mass object of the YSO candidates with optical data. The 
spectral type given in Reipurth & Pettersson (1993) is M2:, so that is the model we have used in Figure [T4l 
It also is one of the stars with a clear B excess in Fig. [131 which can also be noted in the SED itself. Such an 
excess is a characteristic of active accretion. Because it is a previously identified YSO, it appears as grade 
"A+" in Tabled it is a Class II. 

A.5. 073106.5-465454 

073106.5-465454 is the northernmost YSO candidate within the SalOl association. It has a 2MASS 
counterpart, but the SED (Fig. [H)) seems somewhat "disjoint" between the 2MASS and IRAC portions. 
This could be from very high extinction, or intrinsic stellar variability between the epochs of observation 
at the NIR and MIR, or source confusion. It is detected at 24 jim but there is only a limit at 70 /im. We 
categorize it as a "C" -grade YSO candidate, with a "flat" SED class. Additional data are needed to confirm 
or refute the YSO status of this object. 



A.6. 073108.4-470130 

This object, 073108.4-470130, has an SED (Fig. [Ml) quite consistent with other known YSOs, although 
there are no optical data and just a limit at 70 /im. If it is really a young star, it has a significant excess in 
at least 8 and 24 /im, and probably 5.8 /im as well. We categorize it as a very high quality YSO candidate 
("A"), SED Class II. 

A.7. 073109.9-465750 

Like the prior object, 073109.9-465750 has an SED (Fig. [T4| quite consistent with other known YSOs, 
again without optical data. There is a limit at 70 /im, and it is in a high background region, so there is little 
constraint to the SED. Like the prior object (073108.4-470130), if it is really a young star, it has a significant 
excess in at least 2 bands. We categorize it as a very high quality YSO candidate ("A"), SED Class II. 

A.8. 073110. 8-470032=CG-Ha3 

CG-IIa3 is another YSO identified in Reipurth & Pettersson (1993), and has some photometry reported 
there as well. It is detected in all bands considered here, B through 70 /im. In the optical (Fig.fTTl). it is well 
within the clump of YSO candidates above the 30 Myr isochrone. As seen in its SED (Fig. [T^ . the optical 
and NIR data from Reipurth & Pettersson (1993) are quite consistent with the optical and NIR data we 
report here, with some weak evidence for variability in the optical. Reipurth & Pettersson (1993) report a 
type of K7, and that is the model we have used in Fig. [HI Like CG-Ha2, the SED suggests that there may 
be some S-band excess (most likely due to accretion); in Fig. [131 it appears as one of the objects apparently 
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pushed above the main sequence due to reddening. Because it is a previously identified YSO, it appears as 
grade "A+" in Tabled it is a Class II. 

A.9. 073114.6-465842 

This object can be seen at all four IRAC bands and both of the MIPS bands considered here. It also 
has a 2MASS counterpart, but no optical data. The SED (Fig. [T4)) suggests that there is relatively high 
Ay towards this object. Of the new YSO candidates that appear to have some photometric points on the 
photosphere, this object appears to have the highest Ay. If this object is actually a young star, it has an IR 
excess in four bands. We categorize this as another highest-quality (grade "A") YSO candidate, with SED 
Class "flat." Additional data are needed, including optical photometry. 

A.IO. 073114.9-470055 

Another new YSO candidate, 073114.9-470055 has a 2MASS counterpart but no optical data and just 
a limit at 70 iim (Fig. [TS|) . If it is a legitimate young star, it has an IR excess in more than two bands. We 
categorize it as a grade "A" YSO candidate, SED Class II. Additional data are needed, including optical 
photometry. 

A.ll. 073121. 8-465745=CG-Ha4 

CG-IIa4 is another YSO from Reipurth & Pettersson (1993), and has some photometry reported there 
as well as a spectral type of K7-M0. As discussed above, this object has a high-quality J measurement in the 
2MASS point source catalog, but the HKg measurements were flagged as having low-quality photometry. 
We accepted the low-quality measurements with a large uncertainty and added them to the SED in Fig. [151 
Given the the normalization as seen there, there could be a disk excess beginning as early as 3.6 /xm, but with 
a better measurement at the near-IR and thus a better constraint on the location of the photosphere, the disk 
excess could start at a longer wavelength, and we suspect that it probably does. The optical data we report, 
as compared with the optical data from Reipurth & Pettersson (1993), suggest substantial intrinsic source 
variability, a common characteristic of young stars. In the optical color-magnitude diagram in Figure lll[ 
it is apparently the oldest YSO of the set of YSO candidates above the 30 Myr isochrone. This object's 
apparent intrinsic variability could move it around in this diagram; moreover, the uncertain distance to the 
CG4+Sal01 cloud moves this object above or below the 30 Myr isochrone (as discussed above). This is the 
only previously- known YSO in our set that is not detected at 70 fim. It appears in Fig. 1131 as having one of 
the smallest S-band excesses, but the photometry as seen in Fig. [15] suggests that perhaps its placement in 
Fig. [T3| could be improved. Because it is a previously identified YSO, it appears as grade "A-|-" in Table [3] 
it is a Class II. 



A.12. 073136. 6-470013=CG-Ha5 

CG-HaS is reported as a K2-K5 in Reipurth & Pettersson (1993). It is detected at all available bands 
we discuss here. It appears in Figure [T3| as having a i?-band excess and subject to high Ay- Because it is a 
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previously identified YSO, it appears as grade "A+" in Table |3l it is a Class II. 

A.13. 073137.4-470021=CG-Ha6 

CG-HaG, located very close to CG-HaS, is reported as a K7 in Reipurth & Pettersson (1993). It is 
detected at all available bands we discuss here, B through 70 ^m. The infrared excess appears to start 
around 8 /im, suggesting a possible inner disk hole; more detailed modelling is needed to be sure. It appears 
in Figure fT3l as having a S-band excess and subject to high Ay- It is another "A+", Class II, in Tabled 

A.14. 073143.8-465818 

This object, 073143.8-465818, is detected at optical through 24 jum. It is a high-quality YSO candidate 
(grade "A"), and SED class "flat." In the optical fFig. [TT|). it is the lowest mass YSO candidate without 
a prior identification. It has a clear B excess in Fig. [T^l It is not detected at 70 /zm. Spectroscopy of this 
object is required to confirm/refute its YSO status and obtain an initial guess at its mass. 

A.15. 073144.1-470008 

073144.1-470008 is another high-quality YSO candidate (grade "A", Class II) , detected in the optical 
through 24 /im. It too has a clear B excess in Fig. [131 It is undetected at 70 /im. 

A.16. 073145.6-465917 

073145.6-465917 is the last of the YSO candidates in this list associated with SalOl. It is a high-quality 
(grade "A") YSO candidate, detected at optical through 24 /im, but with only a limit at 70 /im. Based 
on the approximate SED fit in Fig. [TH the disk excess starts at 8 /jm, though additional modelling (and a 
spectral type) are required to be sure. In Fig. 1131 it appears as having a B-band excess and subject to Ay. 
It has an SED Class of II. 



A.17. 073326. 8-464842=CG-Ha7 

CG-Hq;7 is the only previously known YSO in the CG4 region covered by our maps. Reipurth & 
Pettersson (1993) report a K5 spectral type; their optical data are quite consistent with the optical data we 
report. We detect this object at all available bands reported here, B through 70 /im. In Fig. [TTl it is the 
highest apparent mass YSO candidate, though reddening can strongly influence this placement. The disk 
excess (Fig. [T5|) begins at the longer IRAC bands, suggesting a possible inner disk hole. It does not appear 
to have a i?-band excess in Fig. [T31 It is located far from any nebulosity (Fig. [T7)) . as noted by Reipurth & 
Pettersson (1993). It is an SED Class II. 
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A.18. 073337.0-465455 

Object 073337.0-465455 has counterparts at J through 24 /im; there are no optical data available. Its 
placement in the [3.6] vs. [3.6] — [24] diagram (Fig. [5]) suggests that it is at the edge of our detection limits; 
it appears at [3.6] — [24] ~3. 4, [3.6] ~14. It is undetected at 70 ^m. If it is a legitimate YSO, it has an excess 
in at least 3 bands. It is located just off the northern edge of the globule. We place it in the "Grade B" bin, 
and it is SED Class II. Additional follow-up data are needed. 

A. 19. 073337.6-464246 

This object, 073337.6-464246, is detected at optical through 24 ^m, but not at 70 fj,m. It is the candidate 
seen in Fig. [TT] as very low in the optical color-magnitude diagram, and closest to the IRAC photospheric 
locus of points in Fig. [51 There is a marginal excess seen at 5.8 and 8 ^m, and then a larger excess at 24 
yum; as discussed above, excesses seen only (or primarily) in a single point at 24 yum can be a result of source 
confusion with adjacent sources. This source is also located very far from any nebulosity, near the edges of 
some of our IRAC maps (Fig. [TT)) . We place this marginal candidate as a grade "C" ; it is an SED Class II. 
Additional data are needed. 



A.20. 073406.9-465805 

073406.9-465805 is the only YSO candidate of our set located projected onto the globule/elephant trunk. 
It appears to also be projected onto a bright rim at 8 /J,m (Fig [17]), so it is potentially being revealed now by 
the action of the ionization front. There are no optical or 70 //m counterparts, but there are counterparts at 
J through 24 /im. If it is a young star, there is an IR excess at more than 2 bands. We place it as a Grade 
"A," Class II object, and additional data are needed. 

A.21. 073425.3-465409 

073425.3-465409 is a very interesting source. In the [3.6] vs. [3.6] — [24] diagram (Fig. [9]), it is the 
brightest, reddest object ([3.6]-11.7, [3.6]-[24]-8.2). And, it is located right on the "lip" of the globule 
(Fig. \17\ , in a region where YSOs might be expected to form. None of the other "fingers" of molecular cloud 
have apparent associations with infrared objects. After the automatic bandmerging described in ij^j above, 
this object did not appear to have a 2MASS counterpart. However, in looking at the overall shape and 
brightness of the SED, we suspected that it should have a match in 2MASS. Examining the images by hand, 
there is clearly a source at this location visible in 2MASS and the POSS images. At the POSS bands, it is 
distinctly fuzzy. It is also slightly resolved at J; it has been identified with 2MASX J07342550-4654106. The 
position given for this source in the 2MASS point source catalog is 2.1" away from the position in the IRAC 
catalog, and the position from the 2MASS extended source catalog is 1.98" away from the IRAC position, 
both of which are very large compared to most of the rest of the catalog. If it is also slightly resolved at 
IRAC, this could affect the positional uncertainty. Certainly, the structure of the molecular cloud around 
this source is complex and could also have affected the position in the catalog. In any case, manual inspection 
of the images ensures that the object is really the same in the two catalogs, so the JHKg photometry from 
the extended source catalog was attached to this Spitzer source. The SED (Fig. [TB]) rises steeply at the long 
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wavelengths, but it is unfortunately just off the edge of the 70 /im maps. 

Given this objeet's SED and loeation, we have classified this as a high-quality (grade "A") YSO candi- 
date, with SED Class I. Additional data are required to eonfirm or refute this object's status. 



A.22. 073439.9-465548 

073439.9-465548, the last object in our list, has counterparts at J through only 8 fim, and as such, it 
has a sparse SED. (It is off the edge of both the 24 and 70 /im maps.) It is located just to the East of the end 
of the globule, consistent with it having been relatively recently uncovered. The shape of the IRAC portion 
of the SED is a little different than that for other objects in this set; despite the negative overall slope of the 
IRAC points, none of the IRAC points seem to be photospheric (compare to other SEDs in Fig. [TiHlHl) . and 
there is even a very slight negative curvature (the slope between 3.6 and 4.5 /im is slightly shallower than 
that between 4.5 and 5.8 /im). On the basis of experience looking at many hundreds of SEDs for YSOs and 
contaminants (RebuU et al. 2010, 2011), we have some reservations about the shape of this SED; without 
additional photometric data, it is hard to give this object a high grade. We give this one a "C" grade; it is 
a Class II. 



